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Motivated by the stripe developments in cuprates, we review some analytical results of our 
studies of the charge- and spin density modulations (CDW and SDW) in a weakly coupled one 
dimensional repulsive electron system on a lattice. It is shown that close to half filling, in the 
high temperature regime above the mean field transition temperature, short range repulsions favor 
charge density fluctuations with wave vectors bearing special relations with those of the spin density 
fluctuations. In the low temperature regime, not only the wave vectors, but also the mutual phases 
of the CDW and SDW become coupled due to a quantum interference phenomenon, leading to the 
stripe phase instability in a quasi one-dimensional repulsive electron system. It is shown that away 
from half filling periodic lattice potential causes cooperative condensation of the spin and charge 
superlattices. "Switching off' this potential causes vanishing of the stripe order. The leading spin- 
charge coupling term in the effective Landau functional is derived microscopically. Results of the ID 
renormalization group (parquet) analysis away from half filling are also presented, which indicate 
transient-scale correlations resembling the mean-field pattern. Farther, the self-consistent solution 
for the spin-charge solitonic superstructure in a quasi-one-dimensional electron system is obtained 
in the framework of the Hubbard model as a function of hole doping and temperature. Possible 
relationship with the stripe phase correlations observed in high Tc cuprates is discussed. 
Stripe phases recently observed in doped antiferromagnets (cuprates and nickelates)"'^"'^ attract attention to the 
problem of multi-mode instabilities in the interacting electron systems. Numerical mean-field calculations and phe- 
nomenological considerations'*"^ suggest a universality of the spin-charge mode coupling phenomenon in repulsive 
electronic systems of different dimensionalities. Hence, analytical mean-field solutions of the multi-mode ordering in 
ID system could be revealing with respect to the mechanism of the mode coupling. The mean-field solutions would be 
essentially dimensional-independent and stabilized in three dimensional model by inter-layer interaction. This review 
is aimed at discussing some of such solutions. 

Inelastic neutron scattering experiments in doped cuprates and nickelates reveal pronounced dynamical spin fluctu- 
ations centered at the wave vectors (I/2±£, 1/2) and (1/2, 1/2±£) (in units of 27r/a, a is the lattice constant). These 
are accompanied by dynamical 'charge' (crystal lattice) fluctuations at wave vectors (±2£, 0) and (0, zt2e). The char- 
acteristic wave vectors of spin and charge are harmonically related via the doping dependent parameter e, indicating 
that mode-couplings are important. Similarly, the static stripe phases show the same harmonic relations between 
the spin- and charge ordering wave vectors, thus confirming an idea that spin-charge mode coupling is a necessary 
condition for these phases to exist. Recent experiments demonstrate^'^ that at least in cuprates one has also to deal 
with the coupling to a pairing mode. It is proposed'' that the stripe phase might be a quantum liquid crystal state of 
electrons possessing at the same time a charge-, antiferromagnetic-, and a superconducting condensates. 

Since the stripes have to do with a multi-mode instability, the experimental discovery of static stripes and dynamical 
stripe fluctuations in cuprates and other doped transition metal oxides^"'^ put forward the following theoretical 
problem: a theory characterized by coupled modes must be extracted from a (microscopic) model of an interacting 
electron system. 

Stripe-like instabilities do show up in a variety of theoretical approaches, some of them even preceding the exper- 
imental discoveries*'^. However, all of these approaches*"^ have in common that they deal with the strong coupling 
regime at zero temperature, while they rest entirely on numerical calculations. To succeed in the same direction ana- 
lytically, a microscopic model with weakly interacting electron-quasiparticles has been chosen*"*". In such a model, in 
order to get at a mean-field ordering transition one limits himself to systems exhibiting near to perfect Fermi surface 
nesting characteristics. This actually implies the choice of only one dimensional (ID) systems. 

At first glance, one may argue here that the infrared fixed point structure of ID interacting electron systems is 
well understood**-**. Namely, it is well established that the zero temperature Luttinger liquid exhibits algebraic long 
range order in both the spin- and the charge sectors*^. Nevertheless, the evidence is growing*^ that existing picture is 
not yet complete. Namely, it is very probable that close to half-filling and for relatively weak interactions the ground 
state is nothing else than a ID stripe phase, with a long range order which is destroyed by a marginal quantum 
fluctuations of the order parameter. It is actually well established, that on the classical level the ground state in 
this regime is a stripe phase, since the mean-field analysis of Schulz^ for 2D stripes rests on the nesting features 
which relate these results with (in depth) a ID effect. Quantum mechanics is involved, as usual, in admixing of the 
Goldstone modes of the classical state in the ground state, causing in turn the algebraic long range order. 

This review consists of the four main parts. The first two parts deal with just the two "main" charge- and spin 
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harmonics coupling, while the third and the fourth parts allow for all the rest "obertones" selfconsistcntly, thus 
introducing solitonic superstructures of the electronic charge and spin. Before coming to the body of the derivations, 
let us make a short outline of the four parts in a more detailed fashion. 

The novelty of the approach used in*'^ and described below lies in that instead of the conventional expansion in 
terms of a (small) coupling constant between electrons, one expands perturbatively in the coupling strengths between 
the (li£fc;rent spin- and charge modes. For this purpose the classical nature of the dynamics of the collective modes is 
used, which has a range of validity in a temperature interval around Tmf - the mean- field transition temperature into 
the ordered state. Thus, the interaction strength between the electrons translates into strength of mode couplings. For 
relatively weak interactions and "high" temperatures one expects on physical grounds a hierarchy of mode couplings. 
The dominant coupling is found between the fundamental spin- and charge modes, just as the £ spin- and 2e charge 
modes of the cuprates. Higher order couplings are responsible for the "solitonic appearance" of the stripe ground 
state, discussed in the last part of the paper. 

Then, considering the coupled spin- and charge fluctuations above Tmf, it is possible to integrate out spin (SDW) 
fluctuations in the presence of the relevant charge (CDW) fluctuations^'^. An effective free energy of the system does 
posses deep local minima at finite CDW amplitudes with the different wave vectors, but the global minimum still 
corresponds to no CDW at all. 

In the second part, results of the paper^° are presented, which directly demonstrate appearance of the global 
minimum of the free energy at finite CDW (and SDW) amplitude below Tmf when the mutual phase of the coupled 
charge- and spin modes is locked for a constructive interference to occur. Hence, a quantum interference mechanism 
of the stripe phase ordering proposed in-'^" is described. In order to carry discussion somewhat beyond the mean-field 
approach, a new modification of the ID ("parquet") renormalization group technique is described, which was first 
applied in^° to the case of coupled charge and spin fluctuations incommensurate with the crystal lattice away from half 
filling. As one would expect on the basis of the well known results for ID systems mentioned above^^-^^, there is no 
four-leg vertex divergence found away from commensurability point in the CDW or SDW channels alone. Nevertheless, 
introduction of the infinitesimal vertices with the new "Umklapp wave vector" 2£, brought by the incommensurate 
CDW mode, causes divergence of the four-leg vertices, thus indicating transient-scale correlations resembling the 
mean-field stripe pattern. 

The last two parts of the review contain results published in-'^^, as well as the new ones. These concern the 
self-consistent solutions of the Bogoliubov-de Gennes equations for a repulsive ID Hubbard Hamiltonian written in 
the Hartree-Fock approximation at zero- and at a finite temperature. Using methods of the finite band potential 
theory the coupled solitonic charge- and spin superstructures are described analytically in terms of the Jakobi 
elliptic functions, which, in turn, represent solutions of the non-linear Schrodinger equations. It is shown that far 
enough from half filling of the bare single electron band the solitonic superstructures smoothly evolve into two coupled 
CDW and SDW harmonics, considered in the first two parts of the paper. A decrease of the effective mode coupling 
constant accompanies this transition. 



I. TWO HARMONICS APPROXIMATION: CHARGE FLUCTUATIONS IN THE CLASSICAL LIMIT. 

This section is based on the results of the works*'^. Our interest is in the description of the precursor fluctuations 
in the metallic state of a one dimensional system, at temperatures above the mean-field ordering temperature. For a 
single mode instability in the weak coupling regime this is well described by RPA. Upon the approach of the mean-fleld 
transition, the fluctuations slow down and in the vicinity of the transition the characteristic frequency of the collective 
fluctuations becomes less than temperature. It was demonstrated by several authors (see e.g.^^ and references therein) 
that under these circumstances the time dependence of the collective modes can be omitted and instead of calculating 
the full quantum- mechanical trace a classical average suffices. 

Let us now consider the Hubbard Hamiltonian with J7 > 0, written in the fermionic second-quantized operators 

H = tY. AaC3,^ + uY,i\nl-{Stf), (1) 

<i,j>a i 

using fiiihii = jhf — (Sf )^ where h = + fi| is the fermion density and the ^-component of the fermion spin. 
The interaction term can be decoupled by the Hubbard-Stratanovich transformation, 

exp[-j{h1 - {Stf)] = ^ / dpidmiexp[-^{pf + ml) + ipifn -t- 2mi5f] , (2) 
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introducing two auxiliary fields, describing the collective charge- {p) and spin (to) fluctuations, respectively. It is noted 
that at temperatures larger than the mean-field transition only spin amplitude fluctuations matter and difficulties^ 
related to the apparent violation of the global SU (2) invariance by Eq. (2) are of no importance. Conventionally, one 
proceeds by neglecting the charge sector completely, except for the q = Q static component of the charge mode causing 
a shift of the thermodynamic potential. However, since the coupling of the charge- and spin modes is important at 
low temperatures the question arises how to deal with these couplings in the precursor regime. 

Starting with a uniform state, it is not easy to keep track of these mode couplings perturbatively. We assume that 
both the charge- and spin mode condense at the mean-field transition where the dynamics in both sectors slow down 
and the quasi static approximation can be applied to either the spin- or the charge modes, or both. Since the spin 
modes dominate the instability, they should be integrated out as accurately as possible (using RPA) which leaves the 
charge modes to be taken as the static ones. Neglecting subdominant charge-charge mode couplings, the trace in the 
partition functional Z can be taken over independent charge modes with wave vector 5c and amplitude p{5c)-, 



Pi = P{Sc) C0s{6cXi) 



and the partition sum can be approximated as. 



Z — J dp{Sc)dScZ{p{5c),6c) 



(3) 



(4) 



where Z{p{6c), Sc) is the partition function describing the dynamics in the presence of the static charge density waves: 
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<i,j>,a 
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TO,;(r)5f(r)]dr 
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4E^i 



(5) 

(6) 
(7) 



where pi is given by Eq. (3) and n is the average number of electrons per site (< 1), while the chemical potential has 
been shifted, 



U 

M--n. 



(8) 



Summarizing, Eqs. (3-7) are a good approximation when the following conditions are satisfied simultaneously: (i) 
the characteristic frequencies of the charge fluctuations should be less than fc^T, (ii) charge-charge mode couplings 
can be neglected, (iii) The modes Eq. (3) exhaust the partition sum in so far the collective charge sector is involved. 
Conditions (i) and (iii) are both controlled by the slowing down associated with the mean-field transition. Although 
close to the transition it is expected that the fundamental charge- and spin harmonics dominate the thermodynamics, 
this is not necessarily the case on microscopic scales, and condition (ii) limits the validity of our approach to the weak 
coupling regime. Self-evidently, because the most important (charge-spin) mode coupling is treated explicitely, the 
above is a considerable improvement over conventional single-mode weak coupling theory, allowing us to penetrate 
deeper into the intermediate coupling regime. 



A. Mode coupling by Fermi-surface matching. 

In this section we will further elaborate on the formalism. However, the central outcome of our analysis can already 
be qualitatively understood at this point: pairs of charge- and spin fluctuations with special relationships between 
their momenta are simultaneously enhanced in the approach to the transition. These include the £-spin and 2e charge 
fluctuations reminescent of the stripe fluctuations in cuprates. In addition, there is actually an infinite progression 
of these pairs, although the enhancement factors rapidly decrease upon going to higher orders. The mechanism is 
straightforward^: by repeated Umklapp scatterings against the CDW, a non-2A;ir SDW can 'borrow' the momentum 



3 



mismatch with regard to the 2kF nesting vector from the CDW. Obviously, this can only happen if special relations 
exist between the wave vectors of the SDW and CDW. 

Close to the transition, the spin fluctuation can be considered as static, as the static charge density wave Eq. (3). 
Due to the underlying lattice periodicity a spin-density wave is of the form 

S{x) = Sq cos{5sx) sin( — ) = Sq- [sin(a;fc_) + sin(xfc+)] (9) 

with fc± = Tr/a ± (5s. The two terms in the brackets in Eq. (9) have a straightforward meaning. The first one arises 
due to the back-scattering of the fermion with a momentum change of fc_ , while the second term comes from Umklapp 
scattering with a total momentum change 27r/a — fc_ = k+. On the other hand, assuming linear fermion dispersions 
the Fermi-momentum can be written in the proximity of half-filling as, 

7 n 6 S TT 

= -z 77 ; t: = Xd— (10) 

2a 2 ' 2 2a ^ ' 

where Xd is the (hole) doping concentration (at half-filling, kp — it /{2a) ). 

Wo have introduced intentionally three different momenta: 5, the deviation of kp from its half-filled value, and the 
wavectors associated with the modulation of the Neel state {5s) and with the charge modulation {6c), respectively. 
Only in the absence of the CDW, the SDW wave vector fc_ equals 2kF = i^ / a — 6 = -k / a — 6s- In the presence of the 
CDW, however, a whole variety of new possibilities arises. According to Eq. (7 ), the fermion momentum is conserved 
only modulo m x 6c, m = ±1,±2,..., because of the presence of the periodic potential produced by the CDW. This 
adds new umklapp scattering vectors, which are linear combinations of the vectors 27r/a and 6c- Therefore, 2kF 
may now differ from either fc_ or fc+ by an integer number of CDW wave vectors: fc_ — 2kp = m6c', m = 0, ±1, 
k,+ — 2kF = n6c', n = 0, 1, .... Notice that the summation over n is restricted to positive integers to avoid redundancy 
in the counting of the allowed SDW wave vectors. These relations are equivalent to, 

6 — 6s = m6c; m = 0, ±1, ±2, ... 

6 + 6s = n5c; n= 1,2,... (11) 

It follows immediately that both relations are only fulfilled simultaneously for the following "matching pairs''^ of wave 
vectors 6c and 6s, 

26 



n — m 

(12) 

{m + n)5 

6s = 

n — m 

where \m\ < n •,m = Q, ±1, ±2, ... ; n = 1, 2, .... The distinction between the matching pairs in Eq. (12) and other pairs 
5c and 5a, obeying only one of the two relations in Eq. (11) will become apparent in the next subsection. Notice that 
the symmetry relations expressed in Eqs. (11-12) arc not related to the actual values of t and U in the Hamiltonian 
Eq. (1), as long as the fermion dispersions are linear and the two-mode approximation (one SDW and one CDW) is 
valid. In combination, this limits the present approach to relatively weak couplings. We will elobarate this further in 
the next sections. 

B. Free energy of the SDW-CDW fluctuations 

Let us proceed to derive the fluctuation contribution to the free energy, in order to sec how the relations Eqs. 
(11,12) follow from the thermodynamics of the electron system described by Eqs. (3-7). For this purpose, the spin 
fluctuations in the presence of the CDW, as described by Eqs. (5-7) will be integrated out on the RPA level. We 
proceed by calculating the free energy functional for fixed charge modulations, 

n{p{6c),6c) = -T\nZ{p{6c),5c) , (13) 
with Z{p{6c), 5c) defined in Eq. (5) and the total partition sum Eq. (4) becomes, 

Z = j d5cj dp{5,) exp(-/3f2(p(^e), 5c)) (14) 
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In the single-loop approximation, Zm{p{5c),^c) in Eq. (6) becomes 

ZM6c), Sc) = exp{-/3(0o + nrn)} ; « ^ + ^^fe (15) 

i 

= - ^ m(^, a;„)m(-^, -a;„)x(5, a;„) , (16) 

where i^(/i) is the density of states at the Fermi level, and x(9, '^n) the Fourier component of the magnetic susceptibility 
of the system calculated with the Hamiltonian Eq. (7), i.e. in the presence of the charge density modulation pi, Eq. 
(3). m,(q,ujn) corresponds with the Fourier component of the Hubbard-Stratonovich spin field at momentum q and 
Matsubara frequency Finally, flj, is a background contribution which is independent of pi and m. Substituting 

Eqs. (15,16) into Eqs. (5,6) and carrying out a Gaussian integration over the real and imaginary parts of the Fourier 
components m{q,Un) we arrive at the following expression for the free energy functional Eq.(13), 

n{p{d,), Se) = -^^'^(m)) ^ ^(5^)2 + fi, + T7V ^ In (1 - Uxiq, ^n)) (17) 

Dropping as usually^^ all the terms with a;„ ^ 0, the semi-static part of the free energy functional (per lattice site) 
i7 is found, depending on the CDW with wave vector 5c and amplitude p(5c), 

h{p{6c, (5c) = ^iP^n - uv{fi)) + T J2 - Ux./a+M^c),Sc)) + . (is) 

q=±Ss 

(lb absorbs the background contributions and will be neglected in the remainder. The first term describes the 

potential energy cost of creating a CDW in a system with repulsive; interactions, and the second term describes the 
decrease of the free energy due to the spin density fluctuations with wave vectors Ss relative to the wave vector of the 
antiferromagnet at half-filling. 

The static magnetic susceptibility x in the presence of the periodic CDW potential remains to be calculated. For 
linearized dispersion of electrons near Fermi momentum it takes the form (compare'^''), 

Xn/a+Ss{p{Sc),Sc) / dx' [exp{i{Ss - S){x - x')}g++{i(j„; X, x')g — (ia;„;x',a;) 

+ (x^x')], (19) 

counting 2kF from its commensurate n/a value according to Eq. (10). The Green's function 5++ [g ) is the slowly 

varying part of the full fermion Green's function, G, in the Matsubara representation for the right (left)-movers. For 
example, 

G++{iuin; X, x') = exp[ikF{x - x')]g++{iu)n; x, x') . (20) 
For these Green's functions we find, 

UJn{x — x') iU 



g++{i(Jn; x,x) = — exp 



2t 



^hen ""("-^^^ > 
2t 



(jj (x — ic'l 

g++{iuJn;x,x') = 0; when < (21) 

using the shorthand p{u) = p{5c) cos{6cu) for the charge modulation. The corresponding Green's function for the 
left- movers is derived by changing the sign in front of the t appearing in the argument of the exponent in Eq. (21). 
Substituting Eq. (21) into Eq. (19), and making use of the identity : 

00 ^ 
Eexp{-.T(2n+l).}--j^; z>0 

n=0 J 
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one obtains: 



( 



\ 



oo exp 
dz 



iz{5s -5)- sin {5cz/2) cos {5c{x - z/2)) 



+ 



i: 



exp 



dz 



sinh {irTz/t) 

iz{5s + 6) + sin {5oz/2) cos {5c{x - z/2)) 

Oct 



sinh {nTz/t) 



(22) 



where lower cut-off in the integration is taken at the lattice constant a = 1. Now, using wellknown property of Bessel 
functions of integer order: 



n— +00 



exp (i^ COS 0) = i"' exp {in4>) Jn{z) 



and averaging over the phase (position) of the CDW in Eq. (22) we find: 



T 
Ji 



dz 



[cos {z{6s + 6)) + cos {z{6s - 6))] ^ ( 2Up{Se) . 



smh{irTz/t) 



Jo 



Set 



■ sin {5cz/2) 



(23) 



Using now the well known^^ addition theorem for Bessel functions of integer order (J„(a;)), we obtain the key result: 



dz 



sinh {nTz/t) 



2 f UpjSe) 



[cos {z{6s + 6)) + cos {z{5s - d))] 



_ T r 

Xn/a+5s — ^ 

+ Jl (-^-^) [^"s (^("^^s + + "-^c)) + cos {z{5s -5 + n6c)) + cos {z{d, + 6- n6e)) 

n—1 ^ c / 

+ COS {z{6s - 6 - n6c))]] 



(24) 



Coarse graining of Eq. (24) with respect to the thermal length of electron, It = i^T/t, leads to the simplified expression 
(«): 



X,r/a+5,(p('^c),(5c) = ^ln^ 



(25) 



where 5a,jj is the Kronecker symbol. For a zero amplitude of the CDW (p(5c) = 0), Eq. (25) reduces to the standard 
weak coupling result and the SDW condensation temperature follows directly, 



2t 



TsDW = — exp ( - — 



TT 



47rt 



(26) 



The picture changes drastically when the CDW's are allowed to have a finite amplitude. To see what happens, let 
us introduce the distribution function describing the probability of finding a CDW fluctuation at wave vector 5c-, 



W{5c) = I 



exp 



^lip(6c),6c) 



dp{6c) 



lil 



exp 



^ip{Sc),Sc) 



1 -1 



dp{6c) dSc 



(27) 



which was calculated numerically. The effect of the thermal length, It, is effectively excluded in the simplified 
expression in Eq. (25) (coarse graining procedure) and is present in Eq. (24). This will especially be of importance at 
large wavelengths, Sc < tt/It, diminishing the multiple scatterings of the fermions against the CDW. It is noticed that 
in this way only an upper bound to the disordering length is incorporated. Correlations are expected to be further 
reduced by charge-charge mode couplings, etcetera, neglected in the present analysis. 

In the figures 1-2 the numerical results following from Eq. (27) are shown for some representative choices of 
parameters. Namely, in Fig.l : Xd = 0.02, U/t = 2; T = 1.2Tsdw; T = 2Tsdw- The blue dashed (T = 1.2Tsdw) 
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and red dashed (T = 2Tsdw) curves are obtained using Eq. (25), i.e. when damping effects have been neglected. 
The blue solid (T = 1.2Tsdw) and red solid (T = 2Tsdw) curves in Fig. 1 were calculated using Eq. (24), i.e. with 
the thermal length effect of the fermion included. This makes peaks sensitive to the temperature, so that higher the 
temperature the more substantial becomes smearing of the peaks. Comparing dashed and solid curves we notice that 
finite thermal length causes merging of the peaks for small values of Sc (corresponding to Zc < 1) into one broad 
feature, while these peaks remain well resolved on the dashed curves at both temperatures chosen. For reasons which 
will become clear in a moment, the wave vectors are normalized to the Fermi surface spanning vector S, 

Zc-J 

Zs^j. (28) 

Let us first consider the undamped case in Fig. 1. Our main result becomes immediately obvious: the charge 
fluctuations exhibit a highly organized behavior as function of their wave vectors. The large momentum regime is 
dominated by a peak centered at Z^ = 2, corresponding with CDW wave vector 6c = 25 which is reminescent of the 
stripe-charge fluctuation as seen in cuprates. This fluctuation becomes more signiflcant both for smalUer dopings and 
smaller temperatures. 

In addition, a variety of smaller peaks is found which occur at finite CDW amplitudes and at wave vectors which 
are rational fractions of the Fermi-surface spanning factor 5, i.e. Zc = Sc/S = n/m, with n, m integers. Since the 
fluctuations under Zc = 2 are not completely unanticipated, these finite amplitude 'fractional momentum' charge 
fluc;tiiations should be viewed as a qualitative novelty, unique to the present analysis. Because they show up at 
smaller momenta, they are also more susceptible for smearing effects as the comparison with the two solid/dashed 
curves in Fig. 1 and two curves in Fig. 2 show. At the same time, it might well be that the members of this series with 
the largest ZcS will survive. Fig. 2 is important for understanding of the Scilla and Horribda situation created by 
the first two terms in Eq. (18). While the first term suppresses peaks at too high doping due to increase of potential 
(" Coulomb" ) energy cost of a CDW, the second term produces less sharper peaks at too small doping due to thermal 
length smearing effect at small wave vector vahie, 5c ^ x^, of a CDW. In order to decrease smearing effect at a fixed 
doping one has to go to lower temperatures, i.e. smaller Tsdw, by decreasing U/t ratio. Solid (black)curve in Fig. 2 
is obtained for the set of parameters: Xd — 0.001, U/t = 1; T = 1.2Tsdw- Dotted (red) curve in Fig. 2 is the same 
as the "low temperature" (i.e. at T = 1.2Tsdw) solid (blue) curve in Fig. 1, and is drawn for comparison. We see 
that peak at (5c = <5 {Zc = 1) has become well resolved, and both peaks on the solid curve are much sharper than on 
the dotted one. This means that effect of CDW precursors with discrete wave vectors demonstarted in Fig. 1 becomes 
well pronounced in the weak coupling-low doping concentration limit. 

The attentive reader should already have realized that the above reflects the 'matching pair' counting of Section A. 
Recalling Eq. (12), the first few members of this scries are: {5s = 5 ,5c = 25}; {5s = 5 ,5c = 5}; {5s = 5 ,5c = 25/3}, 
and {5s = 5 ,5c — 5/2}, which corrrespond with the sets of integers : m = 0; n = 1,2, 3,4. The charge wave 
vectors are clearly recognized in the figures (except = 25/3, which is not seen due to not sufficient sensitivity of 
the Kroncckcr symbol simulation routine at Zc = n/3, n = 1, 2). 

It is instructive to consider how these 'matching pair' fluctuations arise in the present calculation. The key is that 
the CDW fluctuations with a 'matching' wave vector lead to a selective enhancement of the spin density fluctuations 
with proper momenta. The latter push the free-energy downwards, causing pronounced local minima in Q{p{5c,5c) 
(Eq. 18) which are sufficiently deep to carry appreciable statistical weight in Eq. (27). 

The dependence of the free-energy on the CDW wave vector 5c enters entirely through the generalized Lindhard 
function, Eqs.(24),(25). The latter depends on this wave vector both through the argument of the Bessel functions, 
and through the limitations imposed on the summation over the higher order Bessel functions by the Kronecker 5's, 
matching the SDW and CDW vectors to the Fermi momentum. On the other hand, the dependence on the amplitude 
p{5c) enters both through the first term in Eq. (18), describing the restoring force tending to keep the charge density 
uniform in this repulsive system, and again through the argument of the Bessel functions. The first term ^ Jg in Eq. 
(25) is non-zero for any CDW wave vector and is responsible for the 5c independent background. The texture in Figs. 
1-2 is caused by the higher order Bessel functions and in order to gain some intuition one has to investigate which 
values of Sg contribute to the sum in Eq. (18) at any fixed value of 5c, as well as how differs the actual structure of the 
expression xSs in Eq. (25) for the values of 5c, which either belong or do not belong to any matching pair {5c, ds}. For 
this purpose the manifold of points {5c, 5s}, which give non-zero contribution into xSs in Eq. (25), is plotted in Fig. 3 
in the coordinates Zc, Zs- At each point belonging to any straight line in Fig. 3 : Z,, ± 1 = n x Zc', n ~ 0, ±1, ±2 , 
an argument of one of the Kronecker's symbols in Eq. (25) becomes zero . Only Bessel functions of the order n < 8 
were retained in the sum in r.h.s. of Eq. (25) in order to obtain plots in Figs. 1-2. The reason for this is that high 
order Bessel functions, J^{U p{5c) / (tSc)) , give negligible contributions to the sum when their argument becomes much 
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less than order n. The value of the argument, z = Up{Sc)/(tSc), is limited from above by the condition of the high 
probability of the related CDW fluctuation: Q{p{Sc),Sc)/T < 1. This condition leads, due to the first term in Eq. 
(18), to the limitation: z < 4exp {—w)/{'jTXd\/w), where w = int/U, and the estimates: T Tsdw, Sc ^ S were used. 
For the range of doping concentrations Xd and coupling strengths w, which were chosen for the calculations, the last 
inequality gives n < 5. The coordinates of the crossing points of any two lines in fig. 3 determine the matching pairs 
of CDW-SDW fiuctuations with wave- vectors {6c = ZcS, 6s — Zs6}. In order to see the "matching - non- matching" 
distinction in the structure of the expression for xs^ in Eq- (18), Fig. 4 is presented. It shows how many terms 
are simultaneously non-zero in Eq. (25) at each point, which belongs to the manifold in Fig. 3 and has coordinates 
{Zc, Zs}. The latter manifold was substantially rarefied along Zc (but not along Zg !) axis in order to provide better 
visibility in Fig. 4. Each solid circle in Fig. 4 indicates that the number of non-zero terms at this particular point 
is two, while open circles and squares indicate that this number is one. Again, only Bessel functions of the order 
n < 4 are accounted for. The solid circles in Fig. 4 coincide with the two line crossing points in Fig. 3. We see that 
"columns" of solid circles in Fig. 4 erect only upon the 6c {Zc) values which are members of the matching pairs. The 
open circles in these columns would be also solid if we would not restrict the order of the Bessel functions in Eq. (25) 
to n < 8. The structure of the sum in Eq. (25) described above gives rise to the peaks seen in the W(Zc) curves in 
Figs. 1-2. Finally, the absolute walue of W(Zc) is obtained by normalizing to imity on the limited interval of variation 
of Zc- This interval is limited from above by the applicability of the quasi classical approximation for the Green's 
functions of electron, i.e. 6c ^ tt/o ^ 2kF, or in the equivalent form: Zc ^ l/a;^. On the other hand, it is seen in 
Figs. 1-2 that the background part of W{Zc) stretches over an interval of Zc, which is at least order of magnitude 
wide, and extrapolates to a non-zero value in the "forbidden region" , Zc ^ 1/xd- This makes normalization procedure 
somewhat uncertain. 



C. Discussion: 6c = 2ds, but ds = Xd/2 

Here the main results presented in Section I are summarized and discussed. First of all, we notice that in terms 
of reciprocal units, 2n/a, the wave vectors of the symmetry coupled SDW and CDW fluctuations, 'matching pairs', 
chosen at the end of the previous subsection, could be expressed as: 1/2 ± e , ±2e for the first pair; 1/2 ± e , ±e 
for the second pair; 1/2 ± 3£ , ±2£ for the third pair, and 1/2 ± e/3 , ±2e/3 for the fourth matching pair of the 
wave-vectors {6s, 6c}. Here we have according to Eq. (12): e = a6/{2TT) — Xd/2. The whole picture of the CDW 
fluctuation loses its sense if the period of the wave, 2tt/6c, becomes greater than the thermal length. It ^ at/{TrT) 
(compare^*'^°). Allowing for the fact that the most important values are 6c ^ 6s 6 we find that our results are 
limited to the interval of doping concentrations not too close to half-filling: 

T Tsdw ( 47rA 

a=.»-^>-^~exp(^-— j. (29) 

On the other hand, in order a nonperturbativc nature of CDW potential scattering would come into power and prowide 
discrete features scxai in Figs. 1-2, the potential energy cost of a CDW formation with finite amplitude should not be 
too high. Namely, the argument of the Bessel functions in Eqs. (24), (25) should be allowed to be of order unity, i.e. 
Up{6c) > t6c- Substituting then the resulting value of p: p{6c) > t6c/U, into the potential energy term in Eq. (18) 
and requiring that it would be less than temperature T ^ Tsdw one arrives at the upper bound for Xd' 

Xd<^^jexpl-—\. (30) 

Combining Eqs. (29) and (30) we find an interval of 'allowed' values of the doping concentration Xd, where the 
'matching pairs' effect could be expected: 

f 47rA fu ( 2TTt\ 

exp -— — < < W — exp -— — . (31) 



, U J " V i V . 

As the Xd interval above is defined purely via U/t ratio, i.e. coupling strength parameter, we see that in fact Eq. 
(31) defines a region in the phase diagram in the coordinates Xd, U/t, inside which the effect predicted in this paper 
should take place. 

According to Figs. 1-2 the most probable CDW fluctuation at the temperatures T close to Tsdw has wave- vector 
6c = 26. Also, according to Fig. 3, the lowest order SDW fluctuation coupled to this CDW has wave-vector 6s = 6. 
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For the system of a 2-dimensional array of weakly coupled chains (parallel to the x-axis) this (matching) pair of wave- 
vectors: l/2±£ and ±2e, translates into a 2-D form: the SDW wave vector is now (l/2±£, 1/2) , and the CDW wave 
vector is (±2e,0), which coincides with the stripe-theory prediction in a moderately strong coupling limit, U/At > 1, 
for the square lattice, This coincidence is nevertheless quite limited, in the sense that here we have found the 
most probable configuration of the coupled spin-charge fluctuations above Tsdw (see also^), while theory^^ describes 
the ground state spin-charge configuration. On the other hand, when translated into 2-D form, our present result 
is in accord with the inelastic neutron scattering data in underdoped cuprates^ at the temperatures T greater than 
superconducting transition temperature Tc. Namely, the most intensive magnetic neutron scattering is measured at 
2-D wave-vector (1/2 ± e, 1/2), while scattering of neutrons by the crystal lattice potential (which indirectly reflects 
the presense of a CDW) has occured at (±2e, 0). Nevertheless, as in the moderately strong coupling stripe theory^^, 
our prediction for the doping dependence of e: e = Xci/2, differs by a factor of two relative to the experimental result 
in cuprates : e = x^. 



II. STRIPE PHASE ORDERING AS A QUANTUM INTERFERENCE PHENOMENON. 

The main content of this section is based on the work^°. It is well established, that in the Galilean invariant ID 

system one finds on the mean-field level a single mode Peierls wave^"'^ with wave vector 2kp at arbitrary electron 
density {kp is the Fermi momentum). It is generally assumed that the same holds true on a ID lattice away from 
points of lattice commensuration since then an umklapp scattering of electrons by inverse lattice wave- vector 27r/a 
is ineffective^^'^^. As a result, Peierls state with incommensurate periodic spin density structure (SDW) is predicted 
in repulsive case both on the mean- field as on the single loop renormalization group ("parquet'') level. In fact this 
assertion is incorrect. Here we demonstrate that in ID even in a weak coupling limit, U/t 0, the spin-charge modes 
coupling is relevant at small doping Xd, i.e. close to half- filling of electron band {1 — Xd electron per lattice site). The 
mechanism is genuinely novel^*^ and is related to interference phenomenon involving multiple scatterings of electrons 
by inverse lattice wave- vector and by wave- vector of self-consistcntly generated charge density mode (CDW). Essential 
role of the umklapp scattering makes this phenomenon strongly doping-dependent as distinguished from welknown 
phenomenon of related 2A:f-SDW and 4fcF-CDW fluctuations in a Luttinger liquid^^. 

Recently the Landau free-energy functional was phenomcnologicalUy introduced^'', which contains a leading order 
mode coupling term between the fundamental Fourier components of the spin {Sq) and charge (pk) order parameters 

^ = '^i"p\pk\'^ + ^rs\Sq\'^ + \i[{Sq- Sq)pl + c.c.]+ quartic terms (32) 

with condition imposed: k = 2q. Here the third term is relevant signature of the higher harmonics coupling as well. 
Mean-field calculations perfomed for 2D t-U Hubbard models^'^ show that spin-charge modes coupling is important 
at intermediate to strong couplings range of U/t. What happens in ID case, especially at weak coupling? In fact, in 
this case the mode couplings are important at arbitrary small U/t (at the least if doping is small but enough to drive 
the system away from commensuration), and to make this evident a generalization of the usual mean-field strategy is 
needed, which will be discussed below. 

It is wellknown^^, that formation of SDW with wave vector Q- = 2kF, which connects the opposite Fermi points, 
see Fig. 5, may be regarded as Bose-condensation of electron-hole pairs, Cfc^o-c|._|_Q_ o-l*-^)' with the binding energy 
determined by the Q_ -scattering amplitude (here \0) denotes unperturbed vacuum state of the Fermi system). In the 
presence of the lattice potential, the scattering of electrons is split into Q^- and (5+ = 27r/a — 2fci?-scattcring processes, 
the latter involves Umklapp process. In the commensurate case, i.e. at half-filling, a constructive interference between 
the two processes takes place, see Fig. 5a. In the incommensurate case, i.e. at finite doping: Xd — 2-K/a — Akp, 
the Umklapp based (5+-scattering channel brings electron away from the Fermi point region, thus making minor 
contribution to the electron-hole coupling, see Fig. 5b. Electron scattering by a self- consistently generated CDW with 
small wave-vector Qcdw = Q+ — Q- = Xd may (partially) restore constructive interference between Q^- and Q+- 
scattering processes at some "proper" phase shift <p between SDW and CDW. Hence, free energy of thus formed "stripe 
phase", i.e. SDW and CDW with the fixed phase shift^'^'^'^^ might be lower than that of a single incommensurate 
SDW state^". Finally, electron scattering due to a Akp-CDW potential, unlike due to the long-wavelength — Q_- 
CDW potential mentioned above, would interfere with the SDW-induced 2A;F-scattering process at any doping Xd, 
also in the absence of the lattice potential, i.e. in the "empty lattice" case, see Fig. 5c. 

In the presence of the 2fci?-SDW and — (5--CDW condensates, m{x) and p{x), the single-particle eigcnstates of 
the t-U Hubbard Hamiltonian in the Hartree-Fock approximation can be determined from the Bogoliubov-de Gennes 
equations derived in^: 
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du± U , . U , . ^ . . 

^ ~dx' '2'^^^^'^^ ~ ^"^W^=F = (33) 

Here left- and right-movers representation (for momenta close to undoped Fermi "surface" points ±7r/2a) is used for 
the quasi-particle wave function: = M+(a;) exp {ix/A) + au^{x) exp {—ix/A)] and the wave- vectors arc expressed 

in units of 27r/a. The spin density is: m{xi)/2 = (rno/2) cos (Q-X, — ^) = (mo/2) cos (ea;^ + </>) cos (a;i/2), where 
2e = — Q_ = 2-K/a — Akp = x^- Hence, SDW modulation is: rh{x) = {mo/2) cos {ex + (p), and charge density 
varies as: p{x) = po cos2ex. Thus, only slowly varying fimctions u±{x), m{x) and p{x) arc involved in Eq.(33). 

Now, exploit a mathematical trick used in^^ for the problem of field-induced spin density wave (FISDW) in quasi- 
ID (TMTSF)2X compounds. Unlike in^^, here we consider CDW potential instead of magnetic field induced vector- 
potential encountered in FISDW case. Namely, we express wave functions u±{x) in the Bloch-wave basis of the 
periodic CDW potential U p{x) = Upo cos{2ex): 

u±{x) = Cj. exp i^ikx ^ ill / {2v) J po cos{2ex')dx'^ 
= exp {ikx) exp {^iU po sin {2ex) / (4w£)) 

= c1exp{ikx) I Jo(z) + (e"'2- + (-l)»e-'"2-)) j„(^^)| (34) 

where v = 2t is the Fermi velocity of electrons, and z = U poj {Ave). Here expansion of exp \i sinx} in Bcsscl functions 
J„ of integer order n was used. Since J„>i(z) oc when ^ ^ 0, we retain only zero- and first order Bessel functions 
in the last line in Eq. (34) when Upo/{2ve) < 1. After substitution of u±{x) from Eq. (34) into Eq. (33) one finds an 
algebraic system of linear homogeneous equations for the coefficients Cj.. Corresponding determinant equation defines 
the single-particle spectrum^^: 



Eu = ~±^{k+'-) t;2 + A2 (35) 

, A Unio fUpo 

where: A = — - — f — — 
4 \2ve 

f{z) = Jl{z) -2cos{2<f>)Mz)Mz) + J!{z) . (36) 

In the "electron doping" case the sign in front of e in Eq.(35) and of cos2(^ in Eq. (36) should be changed. The 
physical implication of Eqs. (35), (36) is remarkable. Namely, in the absence of CDW z = 0; hence, /(O) = 1 and the 
gap A equals Unio/A (at 2 ^ the Bessel functions behave as: Jo{z) ~ 1 — z^/4:, and Ji{z) « z/2 + 0{z^)). But, in 
the presence of the CDW, A is enhanced when f{z) > 1. The latter condition fixes sign of cos{2(f>) in Eq. (36). The 
form of Eq.(36) manifests quantum interference between scattering amplitudes of electron in the combined periodic 
potentials of (5±-SDW and matching (Q^ — (5_)-CDW. Indeed, r.h.s. of Eq. (36) for /^(z) could be rewritten as: 
\Jo{z) + Ji(2;)e*('^+'^)|^. This is nothing but interference intensity between amplitudes of electron (hole) scattered 
"back" and "forward" by vectors Q- and respectively, with a phase-shift (j> + n. Due to mismatch between the 
scattering wave- vectors, — ^ 0, the interference vanishes in the absence of the matching CDW potential, since 
then: Ji{po = 0) = and Jo(po = 0) = 1. A simple form of solution (35) is valid in the weak coupling limit, U <^t, 
not too close to half filling, i.e. when Xd A/t. 

Free energy of the system (per unit of length), CI, at finite temperature T{= (3~^), follows from Eq.(35) and the 
Hartree-Fock form^ of the Hubbard Hamiltonian: 

O = (C//8) {ml/2 + pI)- {AT/ttv) In [2 cosh {(3E{0/'^)] d4 (37) 

Jo 



where E{^) = VF+A^, and Eb{^ t) is the upper cutoff of the electron energy. Expansion of Eq.(37), in powers 
of small CDW and SDW amplitudes yields: 

^ U , o ,^ A2 / / '^Eb \ 0.053A2 

0»^K/2 + ,J)--(l„(^)-^^ 

^^"•^f,-i,„(^)U,.VW+ ... (38) 



nv \1.76r 

where Ai = —{U^/ {Snv'^Xd)} In [2Ei,/ {1.76T)], and TOq and po are amplitudes of SDW and CDW harmonics with wave- 
vectors e and 2e respectively. We see that Eq. (38) just recovers phenomenological Landau- Ginzburg functional in Eq. 
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(32), considered in^^. Thus, our derivation reveals the quantum interference nature of the SDW-CDW couphng term 
in Eq. (32), with couphng constant Ai fohowing from the microscopic theory, Eq. (38). Notice, that |Ai| increases 
when doping decreases. 

Here we merely list the main results, obtained by minimizing free energy Eq.(37) with respect to SDW and CDW 
amplitudes nio and po- 

i) Coming from the high temperature limit, A = 0, the SDW-CDW superlattice condenses with cos (2(/)) = — 1 or 
cos {2(j)) = 1 depending on the sign of Xd- Thus, the nodes of the spin density coincide with the minima (maxima) of 
the charge density p{x) in the case of the hole (electron) doping, in accord with the stripe phase topology considered 
in the strong coupling limit 

ii) Dimensionless mode couphng strength, Upo/i^TrtXd), in the effective theory Eqs. (37) and (38), grows up to ~ 1 
at small Xd, below Xo ~ \/t/U exp {—2Trt/U). Formal divergence of Ai in Eq. (38) at — > signals that higher order 
harmonics have to be considered as stabilizing solitonic lattice (comparc^^). While transition to solitonic regime is 
governed by parameter Upo/{2TrtXd), the bare coupling constant, U/t<^ 1, remains small. Simultaneously, transition 
temperature, Tc, monotonically increases from the lowest value Tsdw = 2(7/7r)f exp (— 27rt/?7) at \xd\ » Xg, to the 
highest value Tm = 2{j/'!T)texp{—2TTt/{Uf^)) at small doping, \xd\ < Xq. Here 7 = 1.78, and the maximum value of 
the function p{z) in Eq.(36), = P{zm) ~ 1-5, is reached a.t z = Zm ~ 0.83^*. The increase of Tc is accompanied 
by a substantial increase of the SDW amplitude at zero temperature, see Fig. 6. 

iii) The character of the phase transition changes at Xq from the first order {xd < Xo) to the second order 
{xd > Xo), Fig. 7. The jumps of the CDW and SDW amplitudes at the first order transition temperature are: 
™o ~ XdZ„iTjn\'27rtU /{U fmY' and po ~ 2'KXdZmt/U . Hence, in the Xd < Xo region the mode coupling strength is: 
Upo/{27rtXd) » 0.83 (i.e. not <C 1). Therefore, our results in this region, based on the neglect of the higher order 
SDW/CDW harmonics, might be considered as qualitative rather than quantitative. In the H-nd order phase transi- 
tion region, \xd\ ^ Xo, the order parameters close to Tc behave as: Po ~ 3TTgj-,y^/{xdtU), and nio ~ 12^Tsdw /U; 
in qualitative agreement with^'' (here r = 1 — T/Tsdw)- 

The phase transitions described in ii), iii) above belong to the "spin-charge coupling driven" and "spin driven" 
kinds, in the terminology introduced in^^. 

A. STRIPE PHASE AND "PARQUET" ID RENORMALIZATION GROUP TECHNIQUE 

An important issue for the (quasi) ID systems is the influence of fluctuations. We study it within a single-loop 
renormalization group (RG) scheme, so-called "parquet" approximation^^, which we adjust for the case of the two 
order parameters (SDW/CDW) coupled already on the mean- field level. Conventionally, "parquet" - RG equations 
describe behavior of the two-electron scattering vertices 71 (Oi 72(0) and 73 (Oi accounting for back-, forward- , 
and umklapp scattering of electrons respectively close to the Fermi "surface" points: ^kp- The RG variable, ^, 
is the logarithm of the infrared cutoff of the energy/momentum transfer. It is involved in the (logarithmically) 
diverging corrections to the vertices, which are initially defined in the Born approximation: = 7j(^ = 0) <C 1. 
Within "parquet" approach only corrections of the highest power in ^ are retained in each order of the perturbation 
expansion in each Qi and then summed to an infinite order. Transition of the electron system to a strong coupling 
regime is signalled by divergences of the vertices 7^ (0 at some finite value (where "parquet" approximation actually 
fails). In the case of the Hubbard Hamiltonian at half filling: gi — U / (47rt), i = 1,2,3 and = 27rt/[/^^. Away from 
half filling the umklapp condition for two-electron scattering: pi -|- p2 = Ps + P4 ± 27r/a, can not be fulfilled when all 
the quasi-momenta of electrons (before and after scattering) are close to the Fermi surface. In conventional "parquet" 
theory^^ a strong coupling transition does not occur in repulsive ID system away from half filling. The reason is that 
e.g. in the (hole) doped case the deficiency of momentum transfer: 27r/a — Akp = 2e = Xd ^ provides a "natural" 
momentum cutoff, such that at ^ > = \n.{l/xd) the growth of |7i(0l terminates. In order to probe the system for 
a stripe phase instability in this case, we modify "parquet" treatment by adding infinitesimal probe vertices 7i(0) 
which acquire "starting" values 7i(^d) at ^ = S,d- The vertices 71,2(0 describe "umklapp" scattering with the wave 
vector 2e = Xd ^kp, brought by the incommensurate CDW component of the (anticipated) stripe superlattice; 
while 73(0 is due to combined, (commensurate) lattice- and the CDW umklapp: 27r/a — 2£ = Akp- These vertices, in 
combination with the commensurate (bare) umklapp vertex 53, restore possibility of umklapp away from half-filling. 
Thus, "enriched" RG- "parquet" equations in the interval < ^ < 00 become: 

73(0 = -273(074(0; 74(0 - -4i?e(73(073(0) 

74(0 = -273(073(0; 73(0 = -273(074(0 (39) 
where 74(0 = 7i(0 — 272(0' ^'^d same relation is valid between 74 and 71,2- Diverging solutions of Eqs. (39), see 
Fig. 8, can be expressed in the analytical form in the interval ^ > ^d '■ 73(0 = B cosh (C In | tanhZ)(^ ~ ^o)| + <?^o); 
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74(0 = ±y2(7|-S2); 74(0 = {D/2) coth2D(e - 0); 73(e) = i?C/(\/2 sinh 2D(e - 0)), where ~ Co + 0.5(0 - 
Cd) In (2/|7i(ed)|), and all the constants are determined from the boundary conditions for 7^ and 7^ at ^ = and ^ ~ O 
respectively. 

Notice, that position of the stripe-phase ordering point, 0) shifts logarithmically to infinity when starting values of 
the probe vertices 7i(0) tend to zero. This means that long-range incommensurate order is absent in ID, and stripe 
phase order exists at the transient energy(time) /length scales brought by the probe vertices. 

Summarizing, a quantum interference mechanism of the stripe phase ordering in repulsive (quasi) ID electron 
system is proposed^^. Though ID fluctuations smear away mean-field predicted long-range stripe order, "parquet" 
analysis indicates transient scale stripe-phase correlations away from half-filling. The mean-field solution is more 
relevant in the classical spin S ^ 00 limit, and could be stabilized by a coupling to the lattice deformation. 



III. ANALYTICAL STRIPE PHASE SOLUTION FOR THE HUBBARD MODEL: T = 0. 

In this section an exact analytical solution of the Elartrce-Fock problem at T = for a one-dimensional electron 
system at and away from half- filling is described, as it was first found in^^ . This solution provides a unique possibility 
to investigate analytically the structure of the periodic spin-charge solitonic superlattices. It also demonstrates fun- 
damental importance of the higher order commensurability effects, which result in special stability points along the 
axis of concentrations of the doped holes. Though there is no long-range order in the purely one-dimensional system 
due to destructive influence of fluctuations, real cuprates are three-dimensional, and therefore, the long-range order 
survives in the ground state. 

It is well known^° that at low enough temperatures quasi-one-dimensional conductors may undergo a Peierls- or 
spin-Peierls transition and develop accordingly either a long-range charge- or spin-order. In the case of a discrete lattice 
model the commensuration effects become important^^. Exact solutions of related Hartree-Fock problems led to the 
picture of a solitonic lattice appearing away from half-filling of the bare electron band^^'^^. The nature of a soliton was 
determined then by a corresponding order parameter which was either lattice deformation, or density of electronic 
spin. Nevertheless, recently discovered stripe phases in doped antiferromagnets (cuprates and nickelates)^ have 
attracted attention to the problem of coupled spin and charge order pararmeters in the electron systems. Numerical 
mean-field calculations'*'^ suggest a universality of the spin-charge multi-mode coupling phenomenon in repulsive 
electronic systems of different dimensionalities. On the other hand, those calculations are bound to use small clusters 
which often makes them inconclusive. 

Hence, we believe, that one-dimensional mean-field solutions contain universal features of the stripe phase, which 
are stabilized in higher dimensions. We use derived here single-chain analytical solutions as building blocks for the 
stripe phase in quasi two(thrce)-dimcnsional system of parallel chains. In this way we have foimd that short-range 
(nearest neighbour) repulsion between doped holes, in combination with effects of magnetic misfit energy between the 
chains, naturally leads to formation of either "half- filled" or "fully filled" domain walls in the low- or high doping 
limits respectively. In both cases these walls separate neighbouring antiphascd antiferromagnetic domains. 

The Hubbard Hamiltonian with the hopping integral t and on-site repulsion U (> 0) may be written in the form 
already presented in Eq. (1): 

H = tY, + UY.{-X- {St A (40) 

Here an identity: ni-[nii = jfif — (S'?)^ is used. Operators fi = + fii and are fermionic density and spin 
(z-component) operators respectively, and a is spin index. Hamiltonian (40) has convenient form for Hartree-Fock 
decoupling in the presence of the two order parameters, i.e. electron spin- and charge-densities, {S^{x)) and {n{x)) 
respectively. Single-particle eigenstates and eigenvalues of the Hamiltonian Eq.(40) in the Hartree-Fock approximation 
can be determined from the Bogoliubov-de Gennes equations derived in^ (see also a review^^): 

^-i-^&z + cep{x) - am{x)ax^ "^cr = £"^(7 (41) 

where a^^x the Pauli matrices, e = E/2t, a = U/AU the Plank constant is taken as unity, and the length is 
measured in the units of the lattice (chain) period a. In these units momentum and wavevector are dimensionless, 
and velocity and energy posses one and the same dimensionality. The vector ^^r^ = (^'^-i-, ^a-) is defined in terms 
of the right- left-movers ^'ai, which constitute the wave function: 

<i„ [x) = (a;)e''^^/2 + C7*<,_ (a;)e-''^^/2 , (42) 
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where a = ±1 for a spin "f and J, respectively. The Fermi-momentum is ])p = 7rp/2, where in the case of half-fiUing 
the average number of electrons per site equals p = 1, i.e. pp = tt/2. The slowly varying real functions •m{x) and p{x) 
are defined as: 

(n(x)) = p{x), (5'^(x)) = m(x)cos(7rx). (43) 

Note that spin and charge densities have the same coupling constant a. This is not a necessary constraint, our results 
remain valid for a more general case a > /3, where a and fi are spin and charge coupling constant, respectively. 
For a discrete lattice this gives: {S^{i)) = (— l)*m(z). The total energy is equal to^ 



W 
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J2e+ [dx^{m'{x)-p^x)), (44) 



where //, is a chemical potential. Next, we introduce p and p as p(a;) = p + p(a;) and / p{x)dx = 0, and pass to a new 
basis 4*1, ^'2, according to^° (in what follows we drop spin index a) : 

^ ±{x) = eyip{Tia j pda;')*i,2(a;), (45) 
Using this basis we can rewrite Eq.(41) in the form of a single complex order parameter 

-i-^cfz + A(a;)CT+ + A*(a;)CT_ ) * = (e - ap)* (46) 
dx J 

where: A(.t) = — am(.x) exp(2z.a J /5{x)dx), 2(T± = a.j. + ±ay and = {'i>i,'^2)- 

It is easy to find from the finite band potential theory^^ all formal single-soliton solutions of the eigenvalue problem 
(46): 

A(a;) = £0 — ikotanhkox] ko = ^^Aq — eg . (47) 

Here Sq is the energy of the localized level, counted from the chemical potential, and 2Ao is the gap in the energy 
spectrum. Notice, that we consider Aq and £0 (or equivalent pair of variables) as two independent variational 
parameters. The reason is that, according to Eqs. (41) and (46), there are two independent mean fields "hidden" in 
A(a;). Each of them should obey a self- consistency equation. First, we derive equation for m{x), using Eq. (42): 

m{x) = (1/2) ^ {^,+ {x)K_{x) + C.C.) (48) 

<7,E<lJt 

which can also be rewritten in the ^'1.2(2^) representation. Definition given in Eq. (45) leads to a self-consistency 
equation for the variable part of the charge mean-field p{x) = {n{x)) — p : 

p{x) = (*f(ar)*i(ar) + *+(a;)*2(a;)) - p (49) 

The wave functions of the continuum spectrum are as follows: 

±£ T £0 + fc + iko tanh kpx 

*i,2 = _ . e , (p\j) 

2\/L^e{£ - £0) - ko/L 

where £^ = Aq + fc^ , and the upper (lower) sign on the r.h.s. corresponds to the index "1"("2"); L is the length of 
the system (in units of the lattice spacing). 

Simultaneously, the wave functions of the localized state with energy £0 become: 



For correct summation over the energy levels in all the equations above, periodic boundary conditions are imposed 
on the wave functions of the continuum spectrum: '^±{x + L) = '^±{x). Then, quantization condition follows: 
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kL + arctan(fco/fc) = 27rn, where n is an integer. Resulting self-consisting charge and spin components of the single- 
soliton state are obtained: 



p{x)= % r, r^,.-l + -arcsini^, (52) 
2 cosh kqx tt Aq 

m(x) = _ ^0 tanh kox {±2ia0(x) + ix}, (53) 
a 

where 0(a;) = pdx', v is the filling factor of the localized level {y = 0, 1,2); and x is an arbitrary phase. As is 
evident from Eq. (43), at half filling, i.e. in the commensurate case p = 1, the order parameter m{x) must be real. In 
order to fulfil this condition in Eq. (53) (to the lowest order in a for the case = 0,2, and precisely, for v = \') one 
chooses: eq = 0. Eqs. (52), (53) describe structures of the topological spin-charge solitons (kinks), which are either 
spinless with charge ±e (single electron charge) a,t u = 0,2, or chargless with spin 1/2 in the case v = 1. In all the 
cases there arc two antiphased antiferromagnetic domains in the system. The z/ = 1 soliton is the stationary excited 
state of the undoped system. The same is true in the cases v = (j/ = 2), but with one electron removed(added) 
from(to) the system (i.e. still zero doping in the thermodynamic limit). Important is that when doping with holes or 
electrons becomes finite, already the ground state of the system possesses periodic "chain" of the alternating spinless 
solitons and antisolitons of the kind = ot v = 2 (for hole- or electron doping respectively), which are described 
above. This solitonic superlattice, as we assume, is a one-dimensional analogue of the stripe phase observed in lightly 
doped cuprates and nickelates^"^. Existence of the superlattice is related to appearance of the central band around 
e = 0, filled with either spinless "holes" ( = 0) or spinless "electrons" ( v = 2). 

In order to find a structure of the ground state at finite doping nj, = |1 — p| 7^ we calculate the total energy Eq. 
(44) of the electron system using solutions (52) and (53). We obtain for the potential energy, i.e. the second term in 
Eq. (44): 

= ^ - ^ - f r^fco - f V. (54) 

2t 2a a o 2 
where F is defined in Eq. (52). The other part of the total energy in Eq. (44) reads: 

— = (pF£F + Aoln — T ) + w, (55) 

where £% = pjp + A^, and the term w is of the order L°. 

Minimization of the total energy with respect to Aq in the order cx L gives the usual result'^^ : = 2ep exp(— 1/A), 
where A = 2a /n. Since our approximation makes sense when Aq <C £f, we conclude that parameter a = U/At should 
be much less than 1 (weak coupling limit). In this limit we see that part in the potential energy / ap^/2 is of the 
order a and is much less than other terms, which are of the order a~^ or a°. 

The energy of the kink state, Eq. (44), can be expressed similar to the work^^: 

^ = Ao [7 cos 61 -h - sin 61] - ^ A072 sin 6», (56) 

zt TT O 

where 7 = — 26'/7r, and cos^ = eo/Aq, sin^ = fco/Ao. As we have already seen, £0 = for the half- filled case, thus 

leading to 9 = Tr/2. 

In the lowest order expansion in a <C 1 the solitonic structure is described in terms of elliptic functions, see Fig. 9: 

m{x) = {Ao/a)y/qsn{Aox/y^,q), (57) 
p{x) = K{r')Aor [(am(x)/Ao)2 + C] /tt, (58) 

where C = 1 — 2(1 — E{r) / K{r))r~'^ , and su^Aqx / ,Jq, q) is the Jakobi elliptic function with the parameter Q < q <1 
defined by 2K{q)^/ A^ = l/\p— 1|. Here K{r) and E{r) are complete elliptic integrals of the first and second kind 
respectively, and r = 2^/ {q+\), r' — yl^r^. Parameter q varies from g = 1 at p = 1 where sn(Aox) = tanh(Aoa;), 
to g < 1 where sn{Aox / y/q, q) ~ sin(7r|p- l|x). Simultaneously, E{Q)/K{Q) = 1, and E{l)/K{l) = 0. 
In the limiting case of "overdoping" : |p — 1| ^ Aq, in which case 5 <C 1 and K{q) -k /2, one has: 

m{x) « Al{-Ka\p - sin(7r|p - l|x), (59) 
p{x) w A^(7r2|p- l|3)-icos(27r|p- l|a;), (60) 
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in qualitative accord with the approximation of the main harmonics couphng used in-"^*^. Minimization of (56) with 
respect to in the case /? 7^ 1 (it is assumed that |p — 1| ^ Aq/wf, and the ground state structure is described by 
Eqs. (59), (60) leads to the excited state solutions, which are the same as for the Peierls model and are independent of 
a when a <C 1 . The only nontrivial solution is a solitonic excited state with: v = l,eo = Q,Q = Q,W — Wq = ^Aq/tt 
; where Q is the charge of a soliton. This solitonic excitation corresponds to a gradual phase-shift by 7r/2 of the 
argument of sine in the ground state solution in Eq. (59). Thus, we see that the structure of the ground state 
of our model should be similar to the one of the Peierls model^^. Hence, we can conclude that for the finite hole 
density rih = |p — 1| the ground state charged (i/ = 0, 2 ) spin-solitons form a periodic (super)structure with the 
spin period I = 2/\p — 1| equal twice the charge period. It is known that in any exactly integrablc model there is no 
commensurability effects at the commensurate points, : |po ~ 1| = m/n, where m, n are relatively prime integers. 
That is the energy and other system parameters continuously depend on the filling factor p. But in our case the term 
—ap'^/2 in the potential energy (44) violates exact integrability and high order "umklapp" processes lead to a pinning 
of the spin density wave m{x). As a result at any commensurate point we obtain a decrease in the total energy of the 
order^^: 5w oc — aexp{— n const}. 

To summarise, we consider briefly the two-dimensional (2D) case using our ID results by adding weak interchain 
interactions (see also^^ for the charge solitons at half- filling caused by t_L-band effects). In the lowest order approxi- 
mation in the interchain hopping integral t_\_ the interaction energy is 



where J ^ t\/ Aq, ipi is the phase of a spin-density m[x) on the i-th chain, Wc is the Coulomb interaction energy 

between the charged kinks (solitons). We suppose for simplicity that the Coulomb interaction decreases rapidly 
with the distance and take into account only charged kinks on the neighbouring sites of the neighbouring chains, 
Wc = NpQ'^/ea', where Np is the number of pairs of the charged kinks, £ is a dielectric constant, and a' is an 
interchain spacing. In the half-filled system, i.e. in the absence of the charged kinks, the minimum of (61) is achieved 
when (fi — (pj for all i, j, and we have SW = 0. There are two possible ways to create a periodic solitonic structure 
(solitonic superlatticc) in the doped 2D system. In the first case the charged kinks would reside on every chain, 
while in the second case only every even (odd) chain would be doped. Compare the energies of these two different 
configurations. 

In the first case we have (fi = (pj and : 6Wi ~ Wc = NfiQ"^ /2ea' , where Nh is the number of kinks (solitons) and the 
charge Q could be deduced from Eqs. (52) and (58). This is a 2D stripe pattern with filled "domain walls" , i.e. having 
one spinless charged kink per period a' perpendicular to the chains direction. In the second case the Coulomb repulsion 
energy is negligible, but there is an increment in the total energy due to the first term in (61), i.e. due to a magnetic 
order misfit: (fi 7^ (pj. The minimal energy configuration could be achieved now in two ways. One possible pattern 
corresponds to a half- filled "bi-stripe" pattern, seen experimentally in some doped manganites^"^. Namely, the kink 
- anti-kink pairs of the smallest possible size ^ would form on each even(odd) chain in order to keep the phase-shift of 
m{x) with respect to the odd(even) chains equal to 27r. The odd(even) chains remain antiferromagnetically ordered. 
Then the energy is: SW2 ^ NhJ^, where ^ is of the order of the kink width. ^ monotonically increases^'* from its value 
Vf/Aq ~ e^/-^ at p = 1, to vf/A ~ e^/^ tan7r|p - l|/2 in the limit |/0 - 1| » e~^/^ . In the case 6Wi > SW2 we have 
that the alternating kink structure is energetically preferable. Notice, that function dWi monotonically decreases, 
while function SW2 monotonically increases as the function of \p — 1|, therefore the sign in the above inequality may 
change at some sufficient doping concentration of the holes. 

An alternative, half- filled single-stripe pattern may arise when » 2A/7r , i.e. when t± is not too small. In this 
case the minimum of the energy (61) is achieved at ipi = ipj. For this to be true for any couple of the chains, while 
Wc ~ being also fulfilled, the doped holes should again reside, say, on every even chain in the form of spinless 
charged solitons = . But simultaneously, an equal number of the chargeless solitons v = 1 (with spin ±1/2) must 
be formed at all the odd chains in order to maintain the condition (pi = tpj. This configuration will be stable if: 
WgNh < Q'^Nh/{2ea'), where Wg = 2A/n is creation energy of the chargeless kink. Notice that since the charge q 
monotonically decreases from 1 to as function of doping A^^ = Ljl — p| (see Eqs. (52), (60)) the above inequality 
will be not satisfied at high doping densities, and half-filled to filled stripe transition would be expected in qualitative 
accord with experiment 

IV. ANALYTICAL STRIPE PHASE SOLUTION FOR THE HUBBARD MODEL: T > 0. 

Consider thermodynamic properties of the model (1) as a function of a chemical potential ji. The thermodynamic 
potential has the form 




(61) 
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O = -T ^ log(l + e(^-^)/^) + Jdx - J 



(62) 



where p, = p, ~ ap. An analytical treatment is possible at low temperatures T <C Mi ^o, or near a phase transition 
where the gap in the spectrum A <C T. 
We use the solution 



A(x) = Ao\/fcsn(Aox/v^,A:), k = {E+ - E_)/{E++E_), 



(63) 



which is an exact one at the limit a <C 1. The electron spectrum consists of two gaps {E^ < E^ < E^) with the 
dispersion 



dp E"^ - ElE{r)/K{r) 



dE 



E± 



Aq 
2y/k 



(l±fc), r 



k + 1' 



(64) 



where K{k), E{k) are the elliptic integrals of the first and second kind. Using the explicit form for the wave functions^^ 
we rewrite the last term in (62) as 



a 



-2" -2^ 



mdE\ [A4_(A-^)^], 



(65) 



where / is an average value of the periodic function J{x) over a period. After calculations we obtain 

(4-2fc2)£(fc) £;2(fc)- 



A4 - (A2)2 = 1^ 



3 



?,K{k) 



K^{k) 



(66) 



Consider the case of low temperatures T <^ Mi Aq- The thermodynamic potential can be expanded in the parameter 
Ag/T, where Aq is the value of the order parameter at T = fi = 0. Calculating with the density of state (64) we 
obtain after minimization over Aq first terms of expansion of the thermodynamic potential 



2A2 

n = no + — - 



1 

"i 



T 



3/2 



1 



K{k)' K{k) 2{K{k)f 



(67) 



where K « log(4/A;') at A; ^ 1, A:' = Vl — k"^, and 6 = {jl — 2Ao/7r)/(2Ao/7r). The period of superstructure 
(the distance between nearest stripes) is I = 4:K{k)Vk/ Aq. At low chemical potential /U < /ii(T) the minimum of 
the thermodynamic potential is achieved at homogeneous phase (A(x) = const, k = 1). With increasing of fi the 
homogeneous phase becomes unstable, the line of phase transition to the stripe state ( A; ^ 1, I < oo ) is given by 
equation 



ira 
^+24 = 



2 VAo 



3/2 



exp 



Po 



po = 2Ao/-K + ap. 



(68) 



The transition is the phase transition of the first order due to the negative term oc \/K'^ in the thermodynamic 

potential (67). 

The boundary line 12 (/z) between the normal and homogeneous spin density wave state is determined by taking 
fc = 1, Aq — > in self-consistent equations dSl/dk = 0, dSl/dAo = 0. In the lowest order in a this line is given by the 
equation 



log Be 



ilJ_\ 
2 ' 2ttT2J 



0, 



(69) 



where Tq = jAq/tt, 7 is the Eulcr's constant, and is the Euler's digamma-fimction. The phase boundary T^i^) 
between normal and stripe phase is determinated by setting Aq, fc at constant period of superstructure / = 
27r\/fc/Ao. 



log ^ + -Re 



2 



2 27rT3 



1 b-p 



27rT3 



= 0, 



(70) 
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where the wave number of superstructure 2b is the solution of the equation 



At large chemical potential Aq we find T3 oc 

Note that our results are similar to obtained in Refs.^^'^^'^'' for models of one- dimensional superconductors or charge 
density waves. As distinct from previous cases, our model is not exactly solvable, therefore equations (69) - (71) are 
valid in the limit a <C 1. 

The qualitative phase diagram is depicted in Fig. 10. The phase transition between stripe and antiferromagnet state 
is first order, other transition lines are second order. All three lines of phase transitions are intersected at the Lifshitz 
point P. 



= 



(71) 
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VI. FIGURE CAPTIONS 



Fig.l: Probability W{Zc) as function of (normalized) CDW wave- vector Zc = Sc/S; solid blue (dashed blue) lines: 
thermal smearing included (excluded). Parameters: 2t/U = 1, doping concentration: Xd = 0.02. Blue solid/dashed 
lines: T = 1.2Tsdw7 red solid/dashed lines: T — 2Tsdw- 

Fig. 2: Probability W{Zc) as function of CDW wave-vector Zc with the thermal length effect included. Parameters: 
T = 1.2TsDW, "^t/U = 2,Xd = 0.001 - solid (black)line; T = 1.2Tsdw, "^t/U =l,Xd = 0.02 - dotted (red) line. 

Fig. 3: Manifold of points {Zc,Zs}, which give non-zero contribution into X<5^ (p(<5c), <5c) in Eq. (25). The dashed 
vertical grid is guide for the eye. See text for details. 

Fig. 4: Same manifold as in Fig. 3 but rarefied along Z^ axis for better visibility. Solid circles indicate points with the 

number of non-zero terms in Eq. (25) equal to 2; open circles and squares indicate points with this number equal to 
1. The vertical grid is guide for the eye. The "columns" of circles distinguish Z^ - values, which are members of the 
matching pairs. See text for details. 

Fig. 5. Different electron scattering processes in the a) commensurate, b) incommensurate, and c) "empty lattice" 
cases. Lattice spacing is: a = 1, and e = Xd/2. 

Fig. 6. a): SDW (mo) and CDW (po) amplitudes as fimctions of doping at zero temperature, calculated using 
Eq.(37) with At/U = 3.2; b): normahzed stripe phase transition temperature Tj, as function of doping. Regions of 
the I-st and Il-nd order phase transition are separated by vertical dashed line. 

Fig. 7. Calculated temperature dependences of the stripe phase order parameters, nio (curves labeled with "m" ) and 
Po (cmvcs not labeled) for the different doping concentrations x, with At/U = 3.2. Each pair of lines of the same type 
shows TUo and po for each particular value of doping x. 

Fig. 8. Solid lines: numerical solutions of the Eqs. (39) for the vertices 73,4 (curves 3 and 4 respectively) and 73,4 
(curves 3'and 4' respectively), for two different values of ^4: ^1 < ^2- Dashed lines: behavior at half-filling. Initial 
values: 51,2,3 = -94 = 0.1; 73(Ci) = -74(^1) = 0.01. 

Fig. 9. Solitonic spin-charge coupled superstructures for a = U/{At) = 0.3. Only the first harmonics survive at higher 
doping n/( = |1 — p|, in accord with^°. 

Fig. 10. The qualitative phase diagram: the phase transition between stripe and antiferromagnet state is first order, 
other transition lines are second order. All three lines of phase transitions are intersected at the Lifshitz point P. 
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Probability of CDW fluctuations with different wave vectors Zc; solid (dashed) lines: thermal smearing included (excluded) 

Parameters: c = 0.02; U/t = 2; T=1 .2Tsdw (blue); T=2Tsdw (red) 




FIG. 1. 
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Probability of CDW fluctuations with wave vectors Zc; T=1 .2Tsdw 

Dotted line: c=0.02, 2t/U=1, L(therm)/Lc=4; solid line: c=0.001, 2t/U=2, L(therm)/Lc=1 




FIG. 2. 
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FIG. 3. 
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FIG. 8. 
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FIG. 10. The qualitative phase diagram: the phase transition between stripe and antiferromagnet state is first order, other 
transition lines are second order. All three lines of phase transitions are intersected at the Lifshitz point P. 
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